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^ ! Abstract 

Geometric features (including convexity properties) of an exact interior 
gravitational field due to a self-gravitating axisymmetric body of perfect 
^ ' fluid in stationary, rigid rotation are studied. In spite of the seemingly 

00 I non-Newtonian features of the bounding surface for some rotation rates, 

\l . we show, by means of a detailed analysis of the three-dimensional spatial 

' geodesies, that the standard Newtonian convexity properties do hold. A 

central role is played by a family of geodesies that are introduced here, 
and provide a generalization of the Newtonian straight lines parallel to 
the axis of rotation. 



o 

O^' PACS numbers: 04.20.Jb, 04.40.-b 
U 

1 Introduction 

■ In the (thus far, unfulfilled) quest for a realistic exact solution in general rela- 

^ , tivity, representing both the exterior and interior gravitational field generated 

by a self-gravitating axisymmetric mass of perfect fluid in stationary rotation, 
the detailed analysis of the features of whatever partial results we already have 
seem relevant. Specifically, comparison with the known results in the Newtonian 
domain will improve our intuition within the general relativistic regime. 

It is remarkable that there exists a number of treatments based on numerical 
integration of the field equations, or on approximation schemes valid for small 
rotation rates (applied, in particular, for the calculation of the shape of the 
bounding surface of the fiuid configurations, or for the analysis of the meaning of 
centrifugal forces), but, surprisingly, very few exact results, based on the growing 
wealth of interior exact solutions in the literature (both rigidly and differentially 
rotating). In the present paper, we analyze some geometric features of one 
such interior solution, in order to check whether the analog of some Newtonian 
properties hold. Remarkably, they do, in spite of the fact that the analysis of 
the bounding surface p = in Section 3 naively seems to point to the contrary. 
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A more detailed analysis of the three-dimensional geodesies (Sections 4 and 
5) shows that standard Newtonian features indeed hold for the solution under 
consideration, if some of the Newtonian elements are redefined appropriately. 
In particular, we introduce in Section 5 what we believe is the generalization 
of straight lines parallel to the rotation axis in the Newtonian case: Geodesies 
whose points have constant azimuthal angle and intersect the equatorial plane 
orthogonally. 



2 The Kramer solution 

An exact solution of the Einstein field equations, representing the interior grav- 
itational field of a self-gravitating, axially symmetric, rigidly rotating perfect 
fluid, was introduced in Q , and was further analyzed in Q . The metric can be 
written as 

2mds^ = [r;-l-6cosCe"'']dt^ + [4(7?-l)-45cosC(e"''-e"^)]dtd^ (1) 
+ [4(r; - 1) - 4&cos^(e-" + e^-^ - 2e~i)] V _^ ^ 5^^^' ' 

where m and b are positive parameters. The coordinate t is a time coordinate, 
while if is an azimuthal angle. The spacetime possesses the two commuting 
Killing fields dt and d^p . The axis of rotation is characterized by the equation 
r] = 1, and there exists a discrete symmetry ^ — > — ^. The invariant set under 
this symmetry (i.e., points with ^ = 0) will be referred to as the equatorial plane 
in what follows, and the point with coordinates 77 = 1, ^ = as the center of 
the body. The fluid obeys the following barotropic equation of state 

£ + 3p=— , (2) 

where p is the pressure and e the energy density, and kq is a positive con- 
stant. The dependence on 77 and ^ of the pressure and the energy density is the 
following: 

m 

P = — (l + jy-bcos^e-") (3) 
in 

e = —{l-3r] + 3bcos^e-"^). (4) 
2ko 

It is rather remarkable that the pressure is harmonic in the (77, ^) coordinates: 

Pvv + Pec = ■ (5) 

Due to the minimum principle for the corresponding Laplacian, the pressure 
attains its minimum value at the boundary of the domain of definition in the 
(77, ^) plane. This domain is given by the interior of the region bounded by the 
line 77 = 1 and the curve p{ri,£,) = 0. As a matter of fact, the pressure has 
its lowest possible (negative) value at the center. The boundary value p = 
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(the greatest value) defines the boundary of the object. In spite of the pressure 
being negative inside the body, and growing from the center to the boundary, 
the dominant energy condition is satisfied: 

e > 0, bl < £ . (6) 
It is remarkable that the boundary p — has a relatively simple equation 

6cosC-(l + r;)e" = 0. (7) 

This, and the fact that (as will be shown in Section 4) the integration of the rel- 
evant spatial geodesies can be reduced to quadratures, is crucial in our analysis 
of the geometric features of the solution. 

The parameter b is related to the modulus of the vorticity vector at the 
center by means of 

K^^)^ = y£(& + e). (8) 

When the requirement is made that the metric have the appropriate sig- 
nature, as well as the requirement that dt be timelike and 9^ spacelike, the 
following inequalities result [Q: 



m>0, ??>!, 6cosC>0. 1 -r/-|-fecos^e"'' > 0, (r;-f l)(2e''"^ -e^''"^) > 2 . 

(9) 

We shall refer to the intersection of the equatorial plane with the boundary 
p = as the equator of the body, and to the region with ^ > (respectively, 
^ < 0) as the northern (resp., southern) hemisphere. Similarly, the intersection 
of the axis of rotation with the boundary p = having > will be called 
the north pole; the intersection with ^ < will be termed the south pole. The 
object is oblate (in the sense that the polar distance to the center is less than 
the distance from one point in the equator to the center), as calculated in 

As b is bounded away from zero, there is no static limit for this solution. 
One is tempted to interpret this feature in the light of the Newtonian result 
that the pressure cannot have a minimum at the center if il^ < 27rp, where fl 
is the angular velocity of the Newtonian fluid body, and p the mass density. It 
would be of interest to find the corresponding result in general relativity. 

Finally, let us mention that the acceleration and the vorticity vectors are 
parallel at the pole, and orthogonal at the equator, as required for symmetry 
reasons. 



3 Geometry of the bounding surface t =const., 

p = 

The spacetime metric (|l|) induces the following metric on the two-dimensional 
surface t =const., p — 0: 
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2mds2 



[4(?? + l)(2e''-i-e2"-2)_8]d^2 



[62-(l + 77)2e2';](ry2_l) 



(10) 

dff , 



where we have used the equation for the surface p = (|^) in order to express 
the two-dimensional metric as a function of the coordinate r]. It is remarkable 
that (for large enough values of the rotation parameter b) the surface possesses 
a region around the equator with negative Gaussian curvature K (Figs. 1-2). 
Our intuition with two-dimensional surfaces embedded in three-dimensional 
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Figure 2. Gaussian curvature on the two-surface 
t — const., p — 0, in tlie equator, as a function of 
b. 



Figure 1. Gaussian curvature on the two- 
surface t — const., p — 0, from the pole 
{rj — 1) to tile equator, for different values of 
tlie parameter b (b = 6.1,6.3,6.5529,7,7.5) 
corresponding, respectively, to the double- 
dot-dashed, dot-dashed, solid, dashed and 
dotted lines in the figure. (2m — 1) 



Euclidean space would lead us to interpret that region as a concave "waist" 
near the equator. This impression is reinforced by computing the length of the 
parallels (closed curves of constant 7], that circle around the surface and can 
be parametrized by means of the azimuthal angle; they are not geodesies, in 
general). The expression for the length of one such parallel, obtained from (|l0|), 
is the following: 



7r[4(?7-H)(2e 



r,-l 



„27?-2 



It is easily seen that the length presents a local maximum for 



277- 



(2?? + 3)e' 



r,-l 







(11) 



(12) 



As a matter of fact, if the length is plotted as a function of rj (Fig. 3), we see 
that it increases monotonically from 77 = 1 (corresponding to the pole, in which 
case the length vanishes) to a maximum at 77 = 1.1716, obtained by solving 
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Figure 3. Length of the parallels (elosed curves o^ rj — const., parametrized by in the surface 
P = 0.) 



numerically From that point on, the length decreases, until it vanishes 

again for the numerical value r/ = 1.3134. Please notice that the preceding 
values, as well as (12) itself, do not depend on the rotation parameter b. The 
dependence on b, however, shows up in the following: the coordinate rj varies 
from 7] = 1 (intersection of the rotation axis with the surface p = 0) to a 
maximum value (corresponding to the equator), obtained by solving (|^) with 
^ = 0. Accordingly, the maximum value of rj is an increasing function of b. 
When b < 7.0077, the corresponding 77 is such that it falls within the left side 
of the curve in Fig. 3, and we have the "normal" situation, where the length 
of the parallels increases from one pole to the equator. But, if 6 > 7.0077, 
then the maximum length for a parallel occurs at an intermediate latitude, and 
it subsequently decreases towards the equator. In the extreme case b = 8.603 
{rj = 1.3134), the circumferential length at the equator vanishes, which could be 
interpreted as the fission of the body along the equator at the extreme rotation 
rate. Due to positivity requirements in the metric, the values 77 > 1.3134 are 
excluded. 

The Gaussian curvature of the two-surface vanishes precisely at 77 = 1.1716. 
This is not accidental, as the curvature has a factor gcptp^n^ and gipip^n = is 
precisely the condition expressed by (12). Thus, the surface is "flat" at the 
equator for b = 7.0077. If b is increased, then a finite region with K < arises 
symmetrically around the equator, including the equator itself. At the extreme 
value b — 8.603, K becomes singular (minus infinity, see Fig. 2) at the equator. 

If the three-geometry where the two-surface t — const., p = is embed- 
ded were Euclidean, we would find that the three- volume t = const, enclosed 
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by the two-surface would not be convex, and, in particular, certain straight 
lines parallel to the axis of rotation would intersect the boundary p — in. 
more than two points (for sufficiently large values of b) , against the well-known 
Newtonian theorems of Lichtenstein, Q, Q, We shall see below, however, 
that a natural generalization of the mentioned parallel straight lines to the true 
three-dimensional Riemannian geometry does preserve the analog of the classical 
results. 

It should be remarked that the closed curves 77 = const, on the boundary 
surface are not geodesies on the surface, except for the particular case where 
9vv,v ~ ^- This can be readily seen by writing the equations for the geodesies 
in the metric mi 



g^ptpip = const. (13) 
2ffwfy + 9rm,nrf' - gvv,^'^^ = (14) 

(where a dot denotes a derivative with respect to length along the geodesic). 

We thus see that the only parallel circles which are geodesies are the two 
(symmetrically placed with respect to the equator) corresponding to 77 = 1.1716, 
when b > 7.0077. When b = 7.0077, the two parallels coincide with the equator 
(and that is the only case in which the equator is a geodesic). 

Finally, it is easily shown (by numerically computing the Gauss-Bonnet in- 
tegral of K over the surface) that the two-surface t — const., p = has the 
topology of a 2-sphere. 

4 Geodesies of the three-dimensional spatial 
metric 

The spacetime metric (^ reduces for t — const, to the following three-dimensional 
metric: 



2mdso^ ^4\r]~l- 6cosC(e-'' + e^-^ _ 2e"i)] V _^ ^ -r^-—d^^ . (15) 

77 — 1 b cos 4 

From (^5|), we find the following equations for the geodesies in the three- 
space: 



[V-I- 5cos^(e"'' + e^-^ - 2e-^)]ip = const. (16) 
2r/ 77^ . . 



7; — 1 (rj — 1)2 focos^ 



e'-4[l-6cose(-e-'' + e''-")]<^" = (17) 



r^e+|^ + P^C^-46sine(e-"-fe"-^-2e-i)vi^ = 0.(18) 
b cos t, b cos t, b cos^ 4 

In particular, geodesies with ip = are characterized by the two equations 
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77 — 1 (jy — 1)^ 5cos^ ^ ' 

^ + « + = »■ (^o) 

It is a rather remarkable feature of the Kramer solution that the geodesic 
equations ([l^)-(^0|) can be reduced to quadratures. Two clearly different cases 
appear, depending on whether ^ = at all points of the geodesic or not. In the 
former case, the integration reduces to that of the equation 

=0, (21) 



rj-l (?? - 1) 
which yields 

1 



s = q^^-[^^ii^-^^^], (22) 

where g = ±1, depending on the sense in which the geodesic is traversed; rji is 
the initial 77 coordinate, and 77/ the final one, and s is the distance along the 
geodesic. In the case ^ ^ 0, we introduce the new variable w = \Jr\ — 1, in order 
to simplify the equations. By dividing eq. (20) by ^, it can be immediately 
integrated once, giving 

>> 

fce-^" , (23) 



cos^ 

where fc is a positive constant. By substituting (^3|) into (^9|), we get 



If It) = 0, we get the system 



ib - ^we"^' = . (24) 
4&e 



w = ^ (25) 



e2 



cos^ 



ke-' . (26) 



This corresponds to the geodesic along the rotation axis. In the generic case, 
w 0, upon multiplication of ( p4[ ) by w we get 

k ^ 

idw - —e'^''"'ww = , (27) 

whose first integral is 

4^2 + ^e--^ = a (28) 
be 
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with a > a constant of integration. One can now express the relation among 
the coordinate w on the geodesic and the distance s along the geodesic by 



= = qds (29) 



(q — ±1), while the equation for the trajectory is given by 



2eq^ ^ ^ ^ =dw , (30) 



where e = ±1. By using ([l5|), we find a = 2m. To summarize, the relevant 
equations can we written as 

» - (31) 

m 1 



C = 2eg^/--^^-v/^e-- (32) 



bea 



where [3 — ^ . 

Eqs. (|l]) and (||) can be expressed as quadratures: 



^ 2eJ-—^^=dw (33) 



Vcos^ V e 



e 



= iJ^^TP^^^ (34) 

(Notice that q and e are signs, which can be chosen so that the distance s along 
the geodesic increases from the initial value s = at the initial point.) 



5 Convexity properties of the fluid body 

Given the results in Section 3, one could naively expect that the distance from 
points in the surface p = to the axis of rotation would decrease (for large 
enough values of b) as the equator is approached. We shall see that this is 
not the case. In order to do that, let us first identify some general properties 
of geodesies starting from some point in the symmetry axis and reaching a 
point in the boundary p = 0. For definiteness, we shall work in the northern 
hemisphere; due to the symmetry with respect to the equatorial plane, analogous 
considerations hold for the southern hemisphere. 

The first observation is that geodesies from the axis of rotation lie in the 
(w,^) plane, with constant ip. This can be seen from eq. (p^: If the value 
w = {rj ^ 1), characterizing the axis, is substituted, then the left-hand side 
vanishes. Therefore, the constant on the right vanishes. That shows that (p = 
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for such geodesies. The relevant equations for the geodesies are then (p^-(pO|), 
whose integrals are given by (^3|)-(Q). Next, we find the meaning of the constant 
/3: It is easily seen, by using the standard Riemannian formula, that the angle 7 
between the axis and the tangent to the geodesic at the axis (defined such that 
7 = for a geodesic starting at the axis and pointing towards the north pole) 
is related to (3 by the following relation: 

cos7=^. (35) 

For a given point in the boundary p = 0, characterized by w — Wf, it is 
found that a geodesic joining it to the axis has a distance s to the axis given by 



2 rf dw 



m ,ln \J\ — cos^ 76^™ 



(36) 



From (p6D, we see that 
ds [2 



,, sin7C0S7 / :r^dw . (37) 

Vm 'Jo (1 -cos2 7e-"'')i 

As a consequence, s is a decreasing function of 7 in the northern hemisphere. 
The minimum obtains for 7 = f- , which corresponds to a geodesic with constant 
^, whose length (|2|) is 

. (38) 

This, being the minimum of the different distances along different geodesies to 
the axis, we will define as the distance to the axis from the given point (w/,^/) 
in the boundary. Let us now look at how the distance so defined varies when 
we consider different points in the boundary. By using the equation for the 
boundary, 

(u'/ + 2)e'"/ = ^cosC/ , (39) 
and the derivative of Wf with respect to along the boundary. 



dwf ( b 



we get 



ds ds dw f 12 { b 



cosO , (40) 



a<: a \l \ / o 3 I COS^/ , (41) 

at,f owf dt_f \ m \e J 2w^^ + 6wf 

thus showing that the distance from a point in the boundary to the axis increases 
monotonically from the north pole to the equator. 

Another measure of the convexity (or lack thereof) of the boundary p = is 
the behaviour of the distances from the center {w — 0, £, = 0) to points in the 
boundary (Fig. 4). Let us denote by w/(7) the w coordinate of the endpoint of 
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Figure 4. Trajectories of geodesies ip — eonst. from tlie center to the surface p — 0, with different 
values of initial velocities, i.e., varying /3 



a geodesic starting from the center with an angle 7 with respect to the northern 
semiaxis. The distance to the endpoint will be given by 



1 



COS^ 76 



-dw 



(42) 



According to the previous equation, the derivative of s with respect to 7 is 



m ds dwf{'j) 
2d^ " 



1 



tu/(7) 



sm 7 cos 76 

(1 — COS^ 76"'^^) 5 



dw . (43) 



1 — cos^ 76 f ^ ' 

We are interested in the growth properties of s at the equator. In order to 
evaluate the derivative we consider the equation for the boundary (| 

and differentiate it with respect to 7: 



b . 
- sm 

6 



0(7)^ = [2-/(7) + 6.,(7)] e-f^^^'-^ . (44) 



The geodesic with 7 = ■§ corresponds to a geodesic along the equatorial plane, 



with — 0. Hence, C/(f ) — 0- Substituting 7 = -I in 



one gets 

= . (45) 
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But the coefficient in front of the derivative does not vanish; therefore, 

9u>/(7) 







(46) 



By substituting this result in eq. (|4^), evaluated at 7 = ^, we get 

= , (47) 



ds 
97 



thus showing that the distance from the center to the equator is a local ex- 
tremum, if the second derivative of s with respect to 7 does not vanish. Its sign 
will then decide whether the equatorial distance is a maximum (corresponding 
to local convexity at the equator) or a minimum (local concavity). The second 



derivative of s with respect to 7 will be given from eq. (43) by 



2 X 



i7sin7 + cos^ 7W/(7)]e ™/ 



dj 



-sin7COS7e-V(7)(i_cos^e-V(7))-i^!!;£M 

07 



+ (sin^ 7 — cos^ 7) 



-3 sin 7 cos 7 



™/(7) 



(1 — COS^ 76^™^)^ 



-dw 



0/(7) 



(1 — COS^ 76"""^) 2 



■dw 



At 7 = -I , the preceding equation reduces to 



2d^ 



"/(f) 



dw . 



(48) 



(49) 



we first consider the 



In order to evaluate the second derivative ^ g^i'^'' 
equation for the trajectory of a geodesic from the center; from (p3D, 

=dw . (50) 



—== = 2W-COS7 / 
Vcos^ V e Jo 



\/l — COS^ 76 



By differentiating (50) with respect to 7, we find 



aC/(7) 



di ^cose/(7) 



„ /6 dwf(j) 
2W-COS7 — 
V e 07 



-If/ (7) 



sin7 



•\/l — cos^ 76 
(1 — cos^ 7e~"'^)i 



(51) 
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Hence 



(?7 



fb /""'^f^) 
V e 7o 



By differentiating (44) we obtain 



(52) 



5^/(7) 



^7 



-sinC/(7)- 



^{[2«;/(7) + 6«;^(7)]eV(7)} 



dw 



-[2./(7) + 6..,(7)]e-/(^'^^, 



^w/(7) 



and, setting 7 = -^^ in 
■90(7) 
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"/(f) ^^^/(7) 



Finally, from (|54|) and (b2) we get 



(53) 



(54) 



(J^tyy ('J/) 



972 



"'/(f) 



Going back to (^), and substituting ^ '^J'^''') 



972 



(55) 



from (55), the following ex- 



pression for the second derivative of the distance is obtained: 



2 



'"/(f) 



e"" dw 



2{-)'^e~'"f '^^^ /•"'/(f) 



/(f)+3u;/(f)yo 



e-*" dw 



(56) 

But, substituting ^ from (||) (with 0=0) 

we find the following 

inequality: 



2(|)2e-V(f) /-/(f) 2[w/(f) + 2]2e"'/(f) /-/(f) ^ 

e dw — , z — — / e dw 



/(f) + 3«;/(f)yo 



/(f) + 3«;/(f) Jo 



> 



2[u; 2(-) + 2]2e-V(f) /--/(f) 2[u;2(^) + 2]2 



i«/(f) + 3t./(|) Jo 
thus showing that 



5^ 



< . 



(58) 
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It should be stressed that (^) does not depend on b. We conclude that the 
distance from the center presents a local maximum at the equator, thus showing 
that our naive expectations from the analysis in Section 3 were unfounded. 

Let us now consider geodesies joining points symmetrically placed with re- 
spect to the equator, (?«o,Co) and (wq, —^o)i and having (p = 0. For symmetry 
and differentiability reasons, such geodesies must intersect the equatorial plane 
orthogonally, with respect to the metric 

^dw' + ± (i) ■ (59) 

m zm Vo/ cos 4 



The orthogonality condition fixes the parameter (3 in (^3|) and ( |34| ) : 

/3 = ' , (60) 

where Wc is the w coordinate of the intersection of the geodesic with the equa- 
torial plane. We shall now consider the portion of the geodesic starting orthog- 
onally to the equatorial plane from [wc, 0) and ending in (wq, Co), whose length 
will obviously be half that of the complete geodesic starting at (wq,— Co) and 
ending in (wq, Co)- The signs in (|T]) and (^) are fixed by the initial conditions, 
giving 

eq = +l . (61) 

In principle, the considered portion of geodesic in the northern hemisphere 
could have ti; > or w < 0. But the latter does not, in fact, exist, as one would 
have 

13 - e-""' < , (62) 

due to the fact that wq < Wc for w < 0: From (|33|), such a possibility is 
incompatible with the equation for the geodesic trajectory. We conclude that the 
unique geodesic joining {wq, —Co) and (u'o,Co) cuts orthogonally the equatorial 
plane at {wc,0), with Wc < wq, and exhibits the monotonic behaviour w > 
and C > [from (|l]) and (|l|)] in the northern hemisphere. 

Such a geodesic is the natural generalization of a straight line parallel to 
the rotation axis in the Newtonian case: Both can be defined as non-twisting 
{ip — 0) geodesies that intersect the equatorial plane orthogonally. We shall 
now show that the geodesic thus introduced is completely contained in the 
three-dimensional body bounded by the surface t — const., p ~ 0. To this end, 
we consider the pressure p as a function of a point in the geodesic; from (^), 

2—p{s) = w^{s) + 2-- cosC(s)e-'"'("^ . (63) 
m e 

By differentiating ( |63| ) with respect to the distances from the starting point 
(wc, 0) along the portion of the geodesic in the northern hemisphere (we denote 
the derivative by a dot), we find 

2— p= (2u.-f 2-cosCwe~"'')u.-t- -sinCe-^-'C ; (64) 
me e 
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but, due to the inequalities w > and ^ > for the northern portion of the 
geodesic, and the fact that 2^0 + 2^cos^we~"' > and ^sinfe"™ > 0, we 
conclude that 

p>0 . (65) 

As a consequence, the pressure along the geodesic increases from {wc,0) to the 
endpoint {wo,£,o)- Conversely, traversing the geodesic in the opposite sense 
[starting from (wojCo) and heading towards {wc,0)] corresponds to decreasing 
values of p. As the pressure decreases towards the interior in the solution under 
consideration, it is clear that a geodesic of the type just introduced which starts 
at a point in the surface {p = 0) in the northern hemisphere has p < at all 
other points in the northern hemisphere. By symmetry, all other points of the 
geodesic in the southern hemisphere havep < 0, except for the final point, where 
p = 0. Thus, a geodesic of the type considered has only two points of intersection 
with the bounding surface, thus maintaining in the present fluid configuration 
the classical (Newtonian) result for the intersection of straight lines parallel to 
the rotation axis with the boundary p = 0. 

We have found that all the criteria we have considered reproduce the con- 
vexity properties of Newtonian configurations, in spite of the peculiar behaviour 
of the boundary surface, as analyzed in Section 3. It is clear that the standard 
Euclidean relations among convex bodies and their bounding surfaces, Q, 

do not hold for a general Riemannian geometry, dynamically prescribed by 
Einstein's equations. 

6 Conclusions 

We have analyzed some geometric features related to the shape and convexity 
properties of a self-gravitating body of perfect fluid in stationary rotation, as 
given by the exact solution Our analysis of the boundary {t = const., p = 0) 
shows that for some rotation rates there appear features that would be inter- 
preted in the Euclidean case as non-Newtonian. However, by looking at the 
behaviour of geodesies in the three-dimensional fluid within the bounding sur- 
face, the analog of the Newtonian results under consideration is obtained. The 
technique to show them is a detailed analysis of the spatial geodesies within the 
fluid, and specifically the introduction of a family of geodesies which generalize 
the straight lines parallel to the axis of rotation in the Newtonian case. 
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